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Abstract: Many supersymmetric models predict new particles within the reach of 
the next generation of colliders. For an understanding of the model structure and the 
mechanism(s) of symmetry breaking, it is important to know the masses of the new 
particles precisely. In this article the measurement of the mass of the scalar partner 
of the top quark (stop) at an e + e~ collider is studied. A relatively light stop is moti- 
vated by attempts to explain electroweak baryogenesis and can play an important role 
in dark matter relic density. A method is presented which makes use of cross-section 
measurements near the pair-production threshold as well as at higher center-of-mass 
energies. It is shown that this method not only increases the statistical precision, but 
also greatly reduces the systematic uncertainties, which can be important. Numerical 
results are presented, based on a realistic event simulation, for two signal selection 
strategies: using conventional selection cuts, and using an Iterative Discriminant Anal- 
ysis (IDA). Our studies indicate that a precision of Am^ = 0.42 GeV can be achieved, 
representing a major improvement over previous studies. While the analysis of stops is 
particularly challenging due to the possibility of stop hadronization, the general pro- 
cedure could be applied to the mass measurement of other particles as well. We also 
comment on the potential of the IDA to discover a stop quark in this scenario, and we 
revisit the accuracy of the theoretical predictions for the neutralino relic density. 
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1. Introduction 

An experiment at the International Linear Collider (ILC) will be able to make many 
precise measurements from which particle properties, and ultimately, the outlines of a 
particle physics model may be inferred. Due to the high statistical precision expected 
at the ILC, the optimization of the systematic errors is of particular importance. We 
have studied one specific example, namely, the extraction of the mass of an hypothetical 
stop squark from cross-section measurements near threshold. We have devised a method 
which reduces most systematic uncertainties and leads to a potentially very accurate 
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measurement of the stop squark mass. This method, however, is general and could be 
applied to other particles produced in an e + e~ collider. 

The method relies on the comparison of production rates at two different center- 
of-mass energies, and knowledge of how the cross-section varies as a function of \/s 
and the particle mass. In simple terms, one measures the yield at an energy close to 
the pair-production threshold, which will be very sensitive to the particle mass, and 
then at a much higher energy, which has little sensitivity. The ratio of these two yields 
retains sensitivity to the mass, and at the same time is insensitive to many potential 
systematic effects 

We have chosen the case of a light scalar top squark with a mass not much higher 
than the mass of the lightest neutralino since production of this particle was already 
extensively studied in an ILC context [0], [|. It was concluded that a conventional 
approach to the measurement of the stop squark mass culminated in an uncertainty of 
about Amjj = 1.2 GeV [J. The new method improves substantially on this result, 
and for a similar scenario, we conclude that the uncertainty will be Am^ = 0.42 GeV. 

For this analysis, we have performed realistic simulations of the signal and back- 
grounds, and used two techniques to separate the signal from the background. The 
first technique is based on conventional selection cuts, while the second employs an 
improved Iterative Discriminant Analysis (IDA) [|J . Furthermore, the hadronization of 
the stop has been included and we have carefully studied the systematic uncertainties 
arising from this and other sources. 

There are theoretical motivations for studying a light stop squark with a mass close 
to the neutralino mass. Specifically, we evoke a scenario within the Minimal Super- 
symmetric extension of the Standard Model (MSSM) which is able to explain the dark 
matter density of the universe as well as the baryon asymmetry through the mechanism 
of electroweak baryogenesis. The existence of dark matter has been firmly established 
by various observations, most notably by the measurements of the cosmic microwave 
background radiation by the Wilkinson Microwave Anisotropy Probe (WMAP) || and 
the studies conducted by the Sloan Digital Sky Survey 0. The known properties of 
dark matter suggest that it consists of primordial weakly-interacting massive particles. 
Within the context of supersymmetry, the best candidate is the lightest neutralino, 
X\, which is generically the lightest supersymmetric particle, and would be stable if 
.R-parity is conserved. 

Another well-established fact which poses a great puzzle for particle physics is the 
apparent asymmetry between the amount of matter and anti-matter in the universe. 
There are several competing theoretical explanations for the origin of this baryon asym- 
metry. One of these relies on asymmetries generated during the electroweak phase 
transition. The hypothesized mechanism is not viable within the Standard Model 
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(SM), but is possible within the context of supersymmetry. In fact, requiring that 
the correct baryon asymmetry is generated at the electroweak phase transition places 
strong constraints on the parameter space of the MSSM |7|, |], || [K|. In particular, the 
lightest scalar top squark t\ must not be heavy, satisfying the bound < 140 GeV 
with concomitant bounds on the mass of the Higgs boson || [Kj. Furthermore, this 
particle is predominantly of the right-handed chirality state. A small mass difference 
between the stop and the lightest neutralino can help to bring the dark matter relic 
density into the proper range due to co-annihilation between the stop and the neu- 
tralino. For this mechanism to be effective, the typical mass difference is rather small, 

—m^o < 30 GeV [jTTH . The dominant decay mode of the stop is t\ — > cxi, resulting 
in a final state with two soft charm jets and missing energy. Previous studies |2[ 
have shown that clean samples of such events can indeed be isolated at the ILC 

This paper is organized as follows. The next section explains the ratio-of-yields 
method in detail. Section |3| describes the tools and methods used for simulating the 
relevant processes and the detector, as well as two methods for selecting a clear stop 
signal. Section | is devoted to a discussion of the experimental systematics, followed 
by Section || which explores theoretical uncertainties. The last section reports the 
results for this specific channel, and shows the implications for future calculations of 
dark matter relic densities based in particle physics, specifically, supersymmetry. We 
comment briefly on the potential of the IDA method to discover this stop quark at 
\fs = 500 GeV. Conclusions follow. 



2. Method 

One way to measure the stop mass would be through kinematic distributions of its final 
state products. However, jet energies are difficult to measure precisely, especially when 
the jets are not energetic. Furthermore, the radiation of gluons and the hadronization 
of the stop quarks complicate the kinematics in ways that are difficult to predict and 
model accurately. These effects make a precise stop mass measurement from kinematic 
distributions rather difficult 0. 

Alternatively, one can extract the stop mass and mixing angle from measurements 
of the cross-section. For example, it has been shown that using measurements with 
two different beam polarization at one center-of-mass energy, both the stop mass and 
mixing angle can be inferred with good accuracy 0. For light stop quarks with masses 
(9(100 GeV), the typical achievable precision is Am^ ~ 1 GeV. However, this tech- 
nique is limited by substantial systematic uncertainties on the measurement of the 
total cross-section, in particular the modeling of stop hadronization and the resulting 
uncertainties in the selection efficiency 
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We propose a new method which reduces the impact of these systematic uncer- 
tainties, and which we describe in this section in general terms. While our explication 
is based on the case of a light stop, the method could be applied to other particles. 
(See, for example, Ref. |12j] for a discussion of the sensitivity to unknown branching ra- 
tios.) The original presentation of this method concerned Higgs production at a future 
77-collider ]T3|. 

We want to extract the mass (Mx) of a particle from measurements of its produc- 
tion cross sections. In order to obtain the best result, two issues must be considered: 

1. optimization of the energy and luminosity for the minimum statistical error, and 

2. reduction of systematic uncertainties. 

The method described here seeks to address both issues in the best possible way. 

The error on the extracted mass (AMj) relates to the cross-section measurement 
error (Act) through 

It is important to keep in mind that the statistical component in Act depends on a . 

For particles pair-produced mainly in the s-channel, the tree-level cross section 
depends on the mass through the phase space, which usually shows up as factors of 
the velocity of the particle: (3 = yT— (Mx/Eb) 2 where E\, = y/s/2 is the beam energy 
- hence, the maximum energy the given particle can have. For the pair-production 
of scalar particles, a cx /3 3 /s, and for fermions, a oc (3/s. These simple rules can be 
modified by radiative corrections, and by beam energy spread, but the basic picture 
does not change dramatically. We can use this to frame the discussion of the statistical 
error. 

It is instructive to minimize AMx as in Eq. (|2.1|) with a simple Ansatz a = 
<Jof3 3 /s. We imagine that Mx is already known approximately, and we want to select 
the beam energy at which to run the linear collider such that AMx is minimized, 
for a given integrated luminosity C and selection efficiency e. One easily finds that 
AMx = (s 2 / 12(JqM (3) Act. Ignoring systematic uncertainties, one might naively expect 
that Act is proportional to \/N ev , where N ev = aeC is the number of selected events, 
which gives us AMx cx yf/3. This surprising result indicates that zero uncertainty on 
the mass is obtained at the point at which the signal cross-section vanishes. 

The fallacy comes in assuming that Act is proportional to \/N cv , which certainly 
does not apply as N cv — > 0, even in the absence of background. The transition from 
a region in which the cross-section is being measured (Act cx \/N cv ) to a region in 



which an upper limit is being set (N ev ^ 3) is discussed clearly in Ref. |L4]]. One must 
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construct a confidence belt in the (Mx,cr) plane, for a given confidence level — 68% 
would be appropriate for a measurement. In the present case, this belt will depend on 
a(M x ), as well as on eC. When the expected value for N ev becomes too small, there is 
no upper bound on Mx, and Eq. ( |2.1| ) clearly does not apply. In addition, an account 
of background estimates, of experimental uncertainties on e and C, and theoretical 
uncertainties on a(M x ) would require that one does not collect data right at threshold, 
but rather at a point which provides a robust signal somewhat above threshold. 

Equation ( |2.1|) still provides a useful guide in the region above threshold. We 
carried out a Feldman-Cousins confidence-belt construction, and obtained the statis- 
tical uncertainty AMx as a function of the difference y/s/2 — Mx, i-e., the energy 
above threshold. Figure p] shows the result, based on the simple assumed cross-section 
a = a p 3 /s, and approximate values for e and £ corresponding to the analysis de- 
scribed in Section |3.2| . As seen in Fig. [1], the uncertainty on the mass, AMx, really 
does decrease as \fsj2 — > Mx, since the sensitivity of aiMx) to Mx improves more 
rapidly than the relative precision on the cross-section, Aa/a, worsens. 




2 4 6 8 

Vs 12 - M (GeV) 

Figure 1: Statistical uncertainty on the mass, AMx, as a function of the beam energy 
above threshold, \fsj1 — M x . This result is based on a Feldman-Cousins confidence belt con- 
struction, with a simple Ansatz for a{Mx) and approximate values for e and C Backgrounds 
were not taken into account. 

We turn now to a more realistic case. The theoretical cross section as a function 
of y/s is shown in Fig. |2], for two stop masses (122.5 GeV and 123.5 GeV). We include 
QED radiative effects, as described in Section [5]. Following the scenario discussed in 
Ref. [Q, we consider Mx ~ 123 GeV, and beam polarizations of 80% for the electron, 
and —60% for the positron. If we want to use a measurement of the cross section to 
distinguish these two masses, then clearly the threshold region is the most sensitive. 
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This corresponds to maximizing da/dMx, which will minimize AMx according to 
Eq. (|2.1|) . The lower plot on the left side of Fig. |2] zooms in on the threshold region, 
to show how much the cross-section differs for two different hypothetical stop masses, 
and the lower plot on the right shows this difference relative to the cross-section for 
m h = 123 GeV. 

Recall the relation of the cross section to experimental quantities: 

< 2 - 2 > 

In a real analysis, B, e and £ all carry systematic uncertainties, which must be assessed 
and taken into account. An 'optimal' analysis will keep these to a minimum. 

Usually the most difficult component in the systematic error comes from the effi- 
ciency and acceptance. An absolute cross section requires knowledge of the absolute 
efficiency, which, in the case of the ti search described in Ref . , involves charm-tagging 
as well as the hadronization and fragmentation of the t\ and c-quark. While a large 
sample of e + e~ —>■ t\t\ events will allow one to tune Monte Carlo models, and other 
Standard Model processes may provide large samples of c-jets for measuring efficiencies 
for c-tagging, it may be useful to have a method which is relatively insensitive to these 
sources of systematic uncertainties. 

The common step toward reducing systematic uncertainties from the efficiency is 
to work with ratios of cross sections. This also can reduce uncertainties from the 
luminosity measurement, and potentially, from the background and theoretical signal 
cross-section as well. We propose to measure the yield of signal events close to threshold, 
which will be very sensitive to Mx, and compare it to the yield near the peak of the 
excitation curve, which will be insensitive to Mx (see Fig. |2|). We define the observable 



Y(Mx, Vs^) = — — = , - — v -^r- (2.3) 

JVpk — -Dpk tJ t{y/ s pk) Cpk *-pk 

where Nth and B t h are the numbers of selected events and estimated background events 
for y^Sih near threshold, and N p ^, B p ^ are the same quantities for -y/Spk near the peak 
of the excitation curve. Anticipating the results of later sections, we have computed 
the observable Y as a function of , and displayed the result in Fig. |[ 

The slope of the line in Fig. |3] depends on several factors, and one can attempt to 
optimize Y in order to obtain the best measurement of . The sensitivity of Y to 
comes through (Ji(^fs^h), so ^fs^h should be close to 2m^, as discussed above. Mindful 
of large theoretical and growing experimental uncertainties as y/sth 2m^ , we have 
selected y / Sth = 260 GeV, which is 14 GeV above the nominal threshold for a stop 
with = 123 GeV. We find the peak cross-section occurs at v /s p T ~ 370 GeV, but 
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Figure 2: Cross-sections for pair production of the lightest stop squark. Top plot shows 
the full excitation curve as a function of \fs for two slightly different values of mf . The lower 
left-hand plot shows a close-up of the threshold region. The lower right-hand plot shows the 
difference of the two cross-sections relative to their average value. Clearly the largest relative 
difference is seen very close to threshold. 

y/s = 500 GeV would also serve well. Reducing the statistical uncertainty on Y to an 
absolute minimum would require maximizing the integrated luminosity at threshold, 
£ t h> but in reality one would not run the ILC at v /ith = 260 GeV for very long, and in 
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practice C t h = 50 fb 1 is already adequate. We assume £ p k = 200 fbf 



We computed the cross sections with the program Calvin ||15|| , which includes next- 
to-leading (NLO) order supersymmetric QCD corrections, and which was modified for 
this work to include resummed Coulomb corrections near threshold (see Section [|). For 
two common choices of beam polarization, the cross-sections are 

P(e~) = -80%/P(e+) = +60% : a(y/s£) = 17.4 fb a(^s^) = 72 fb, 
P( e -) = +80%/P(e+) = -60% : a(y/s£) = 77 fb a(y/s£) = 276 fb, 



where P < stands for left-handed polarization and P > for right-handed polariza- 
tion. We choose the second set of polarization values since it leads to a much better 
signal-to-background ratio. 

For the computation of the observable Y depicted in Fig. |, we employed the results 
of the "cut-based" analysis described in Section |3.2[ The efficiencies at threshold and 
peak are e t h = 0.34 and e p k = 0.21, and the total background cross-sections are 2.5 fb 
and 10.3 fb, respectively. The strong variation of Y with in Fig. |3] indicates that a 
precise measurement of Y will lead to a precise value for mj . The shaded horizontal 
band corresponds to a 3% uncertainty on Y, resulting in Am^ = 0.2 GeV, which would 
be far better than the result reported in Ref. |]. 

We consider now the impact of systematic uncertainties on the observable Y, and 
Eq. (|2.3|) provides our starting point. For the event selection criteria described in 
Sections |3.2| and |3l| the signal is much bigger than the background, so the main 
experimental uncertainties will come from e. The values for e at threshold and on the 
peak come from Monte Carlo simulations of the signal process. Systematic errors arise 
when these simulations do not match reality perfectly. For example, the calibration 
of the calorimeter energy measurement for real data may be slightly different than is 
simulated, in which case the efficiency for a cut on the total visible energy E v i S as 
estimated from the simulation will be slightly incorrect. One can express the impact 
of this error on the efficiency as e true = e est (l + 5), so that 5 is the relative shift in the 
efficiency. Then the impact on the observable Y is simply 

/I i J! \ vtrue vest 

-true _ vest / 1 + °th \ > ~ > 



Thus, if the systematic uncertainties 5 P ^ and <5 t h are correlated, and if they have the 
same relative impact on e p k and e t h, the net effect on Y will be zero, and there will 
be no error on mt . For some systematic effects, the errors will be correlated, but of a 
different magnitude at the two energies, so that the cancellation | <5 p k — 5 t h\ will not be 
complete. For other systematic effects, the errors will be uncorrelated, in which case 
there is no cancellation. Clearly the analysis should be designed in such a way as to 
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Figure 3: Variation of observable Y with m^, shown by the solid red line. The horizontal 
line gives the expected value of Y when mi = 122.5 GeV, and the shaded band shows a 
variation of Y by 3%. The vertical arrows indicate the corresponding uncertainty on mf . 

take advantage of this cancellation. In practice, this means that the cuts should have 
a similar impact on the signal at both energies. For the present application, there is a 
large degree of cancellation, leading to a greatly reduced systematic uncertainty on the 
observable Y, and hence on m^. The details are given in Section [|. 

We proceed now to a detailed and realistic simulation, and the description of two 
fully-developed event selection methods. 



3. Event Selection and Analysis 

At an e + e~ collider, scalar top quarks would be produced in pairs, and decay to a 
c-quark and the lightest neutralino: 

e+e-^M;->cx°cx°. (3.1) 

The stop quarks live long enough to hadronize before decaying, so the final state signa- 
ture consists of two charm quark jets, missing energy and possibly additional jets due 
to the hadronization process and gluon radiation. 
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In the following sections, the method described in Section |2] will be applied to the 
theoretical parameter point of Ref. ||. The weak-scale MSSM parameters are 

m | 3 = _99 2 GeV, wiq 3 = 4330 GeV, 

Mi = 118.8 GeV, M 2 = 225 GeV, 

A t = -1100 GeV, m A o = 800 GeV, 

The corresponding tree-level masses are: 

m h = 122.5 GeV, m h = 4333 GeV, 
77V) = 107.2 GeV, m*o = 170.8 GeV, m,o = 231.2 GeV, m*o = 297.7 GeV, 

Al A2 A3 A4 

77i-± = 162.7 GeV, m-± = 296.2 GeV, 

Xl A2 

(3.3) 

and the light stop state is almost completely right-chiral, cos#£ = 0.010. As a result of 
the small stop-neutralino mass difference, the stop almost completely decays through 
the loop-induced process into a charm and neutralino, t — > c%°. Due to the loop 
suppression of the decay, the stop is expected to hadronize before decaying. We have 
carried out realistic experimental simulations, and will present the analysis of relevant 
systematic effects. 



Q,U,D,L,Ri, 2 
I/' I 



10 TeV, 
225 GeV, 

0.2, tan/? = 5. 

(3.2) 



3.1 Simulation 



Both the signal and background events are generated with Pythia 6.129 [16]. The 



cross-sections for the signal process were computed with Calvin [15] with some im- 
provements as in Ref. [T7|]. The relevant background processes have been computed 
by adapting the Monte Carlo code used in Ref. |DJ and by Grace 2.0 [[LJ|, with 
cross-checks with CompHep 4.4 [IT5]. The simulation and cross-section calculations 



incorporated beamstrahlung for cold ILC technology as parameterized in the program 
Circe 1.0 ||20|| . Table p] summarizes the predicted signal and background cross-sections. 
To avoid the infrared divergence of the two-photon background process, a cut on the 
minimal transverse momentum is applied, pt > 5 GeV. Backgrounds from supersym- 
metric processes will be discussed below. Table ^] lists the numbers of events generated 
and equivalent luminosity based on the cross-sections in Table [|. 

Hadronization of the final state charm quark and the intermediate stop quark are 
a key issue in this study. The Lund string fragmentation model was used together with 
the Peterson fragmentation function plf . The stop fragmentation is simulated [22] 
by labeling the stop quark as a stable particle in an intermediate step, and switching 
on the stop decay again after stop fragmentation. The modeling of the hadronization 
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Process 


Cross-section [pb] at v 


^ = 260 GeV 


Cross-section [pb] at v 


= 500 GeV 


P(e-)/P(e+) 


0/0 


-80%/+60% 


+80%/-60% 


0/0 


-80%/+60% 


+80%/-60% 


Ux 


0.032 


0.017 


0.077 


0.118 


0.072 


0.276 


w + w~ 


16.9 


48.6 


1.77 


8.6 


24.5 


0.77 


zz 


1.12 


2.28 


0.99 


0.49 


1.02 


0.44 


Wev 


1.73 


3.04 


0.50 


6.14 


10.6 


1.82 


eeZ 


5.1 


6.0 


4.3 


7.5 


8.5 


6.2 


qq,q^t 


49.5 


92.7 


53.1 


13.1 


25.4 


14.9 


ti 


0.0 


0.0 


0.0 


0.55 


1.13 


0.50 


2-photon 

p T > 5 GeV 


786 






936 







Table 1: Cross-sections for the stop signal and Standard Model background processes for 
y/sth = 260 GeV and ^/s p k = 500 GeV and different combinations of beam polarization. The 
signal is given for a right-chiral stop of = 122.5 GeV. Negative polarization values refer 
to left-handed polarization and positive values to right-handed polarization. 



spectrum of the stop is described in Ref. [23]. The dominant lightest stop hadron states 



are mesons composed of a stop and an up or down quark. 



The Simdet detector simulation [24] was used, describing a typical ILC detector. 
The analysis used the N-Tuple tool P5 |, which incorporates jet-finding algorithms. 
In order to reduce the size of the ntuples, several pre-selection cuts were applied, as 
was done for the previous analysis [Q: 

4 < V tracks < 50, p T > 5 GeV, 
I cos0 thrus t| < 0.8, |pi/ptot| < 0.9, 

E vis < 0.75Vs, m inv < 200 GeV. (3.4) 
Most of these cuts have very little impact on the signal efficiency 
3.2 Sequential-Cut Analysis 

Although Standard Model background processes are several orders of magnitude larger 
than the stop signal process, the background contributions can be reduced to an ac- 
ceptable level by suitable selection cuts. This work follows the analysis of Ref. jjj], but 
makes some adjustments to accommodate the stop fragmentation effects, and to take 
advantage of the cancellation of systematic uncertainties as discussed in Section |2|. 

The event selection begins with some basic and common kinematic cuts based on 
global event quantities. The visible energy, E v - m , must be less than 0.3-y/i to ensure 
a large missing-energy signature. It must be greater than O.ly/s to suppress the bulk 
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^ = 260 GeV 


= 500 GeV 




generated 


luminosity (fb~ 




generated 


luminosity (fb 1 ) 


P(e-)/P(e+) 




0/0 +80%/- 


60% 




0/0 


+80%/-60% 




50,000 


1562 


649 


50,000 


423 


181 


w+w~ 


180,000 


11 


102 


210,000 


24 


273 


zz 


30,000 


27 


30 


30,000 


61 


68 


Wev 


210,000 


121 


420 


210,000 


34 


115 


eeZ 


210,000 


41 


49 


210,000 


28 


34 


qq,q^t 


350,000 


7 


6 


350,000 


27 


23 


ti 








180,000 


327 


360 


2-photon 


1.6 x 10 6 


2 


2 


8.5 x 10 6 


9 


9 



Table 2: Numbers of generated events at y / ith = 260 GeV and ^/s p k = 500 GeV, and the 
equivalent luminosities in fb -1 . 



of the two-photon events. Similarly, the number of reconstructed charged tracks must 
indicate real hadronic jets, so we require iV trac k s > 5. In order to suppress Wev and qq 
signals, we place an upper bound V trac k s < 25 at threshold and V trac k s < 20 at peak. 
These cuts on V trac k s remove only a couple percent of the signal. 

We place one more kinematic and one topological cut to further reduce the back- 
grounds. The cuts values are carefully tuned to achieve a low systematic uncertainty 
for the observable Y, as well as a good background rejection. In practice, this means 
aiming to remove approximately the same amount of signal at the two center-of-mass 
energies, rather than achieving the highest signal efficiencies. In particular, the effi- 
ciency at y/Spk is relatively unimportant since we anticipate a large luminosity and 
a large signal cross-section there. The thrust value is useful for eliminating qq and 
two-photon events. As shown in Fig. f|, the thrust distribution for the signal is rather 
different at the two center-of-mass energies, so we require 0.77 < T < 0.97 at ^/si^ 
and 0.55 < T < 0.90 at y/s p ^. Similarly, the event pt, calculated from all energy 
flow objects in the event, is crucial for eliminating the two-photon background. Our 
study indicates that a minimum cut pt > 15 GeV is needed. We tighten this cut to 
Pt > 22 GeV at v /s p T m order to eliminate the same amount of signal events as are 
eliminated at ^/s^. Fig. ^| shows that cutting at px = 22 GeV at y/Spk places the cut 
at almost the same point in the px distribution for both center-of-mass energies. An 
upper cut on pt helps reduce the Wev background, so we require pr < 45 GeV at i/sth 
and pt < 50 GeV at v /s^k, which again reflects our effort to minimize the systematic 
uncertainty. 

One might expect that the signal process ( |3.1|) produces only two jets. However, 
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Figure 4: Signal marginal distributions for the thrust, at y/sth (left) and at -^/Spk (right). 
The arrows indicate the cuts chosen to minimize the systematic uncertainties as well as the 
statistical uncertainty. 




10 20 30 40 50 60 

event p T (GeV) 




event p T (GeV) 



Figure 5: Signal marginal distributions for pr, at y / ith (left) and at ^/s p k (right). The ar- 
rows indicate the cuts chosen to minimize the systematic uncertainties as well as the statistical 
uncertainty. 

additional soft jets can emerge from the stop hadronization process and also from the 
decay of the stop hadron. In order to maintain a high efficiency, and to avoid large 
systematic uncertainties from the modeling of the rate and characteristics of these extra 
jets, events with more than two jets should not be rejected. However, to suppress the 
background processes effectively, extra jets are allowed only when their energy falls 
below a certain cut-off value. To be specific, if there are more than two jets in an 
event, only two of the jets are allowed to have energies above 25 GeV. In this paper, 
we refer to this requirement as the "extra-jet veto." Furthermore, if there are more 
than three jets, the most energetic jet cannot be too energetic - its energy must be 
less than 35 GeV. These cuts are useful against the troublesome Weu backgrounds, 
especially at y/Spk- 

Further substantial improvements of the signal-to-background ratio can be achieved 
by exploiting kinematic and topological correlations between the two c-quark jets. 
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Therefore it is necessary to identify them from the plurality of jets, and for this we 
use charm tagging as realized using a neural network ||26|| . The neural network uses 
information about the vertex position of a jet based on a topological vertex finder, the 
impact parameter probability, the momenta of the associated tracks and the recon- 
structed mass. It has been optimized to single out charm jets with an energy that is 
typical for light stops, while rejecting light quark jets coming from Wev background. 
Each jet in an event is tested with the charm tagger, and a charm flag iy is set (ide- 
ally, f/ = 1 indicates a charm jet). First, if a displaced vertex is reconstructed, the 
jet is tagged positively with F^ = 1. A displaced vertex is found roughly 50% of the 
time for a charm jet, and less than 20% of the time for a light quark jet. If no such 
vertex is reconstructed, then the neural network is employed, which produces a charm 
flag value between zero and one, < F^ < 1. The output of the neutral network is 
shown in Fig. BL for the second of the two charm-tagged jets. 

(c) 

We consider the two jets in the event with the highest values of iy , and require 

(c) (c) 

P c = Ff' x Ff> > 0.6, which is very effective at eliminating events with no charm- 
quark jets while retaining a high efficiency for signal events. In particular, the Wev 
background is reduced by more than half. Fig. ^ compares the quantity P c for signal 
events and Wev background which have passed the kinematic event selection cuts. 
Since half of the Wev events have a genuine charm jet, it is the value of iy for the 
second jet which best distinguishes signal and background. 

A further substantial reduction can be obtained from cuts on the invariant mass of 
the the two best charm-tagged jets - we veto events in which that mass is consistent 
with the PF-boson mass. 

0.08 n 




neural network output 

Figure 6: Illustration of the ability of the neural network to discriminate signal and the 
main background coming from Wev production, for the second of two charm-tagged jets. 
Both distributions are normalized to unit area. 
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Figure 7: Event charm probability, P c = F^ x F^' , comparing signal and the Wev 
background. Our requirement is P c > 0.6. These distributions are normalized to unit area. 
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The event selection cuts are summarized in Table ^ for the two center-of-mass 
energies, A/sth = 260 GeV and y/Spk = 500 GeV. These follow the pre-selection cuts 
listed in Eq. fj). 



Variable 


^ = 260 GeV 


ysT^ = 500 GeV 


number of charged tracks 
visible energy E vis 
event longitudinal momentum 
event transverse momentum p? 
thrust T 

number of jets Nj ets 
extra-jet veto 

charm tagging likelihood P c 
di-jet invariant mass mjj 


5 < iV tracks < 25 
0.1 < E wis /yf~s < 0.3 

biMotl < 0.85 
15 < p T < 45 GeV 
0.77 < T < 0.97 

iV jets > 2 
£ jet < 25 GeV 
P c > 0.6 
ml < 5500 GeV 2 or 
ml > 8000 GeV 2 


5 < AT tracks < 20 
0.1 < E v J^~s < 0.3 

|Pz/Ptot| < 0.85 
22 < p T < 50 GeV 
0.55 < T < 0.90 

iVjets > 2 

E iet < 25 GeV 
P c > 0.6 
m 2 j < 5500 GeV 2 or 
ml > 10000 GeV 2 


signal efficiency 


0.340 


0.212 



Table 3: Selection cuts for y / ith = 260 GeV and -y/s^k = 500 GeV. Also listed are the selec- 
tion efficiencies for right-chiral stop squarks and neutralinos with masses given in Eq. ( |3.3| ). 
See the text for explanations of the extra-jet veto, charm tagging, and the m 2 j cut. 

Our estimates of the numbers of signal and background events surviving the cuts 
listed in Table ^| are summarized in Table If, in a given channel, no simulated 
events remain after applying our cuts, we list an upper limit corresponding to one 
simulated event, and we count this amount in the total background estimate. As 
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V^th 


= 260 GeV 


V^Pk 


- 500 GeV 




C- 


= 50 ib" 1 


C = 


= 200 ftr 1 


P(e-)/P(e+) 


0/0 


f80%/-60% 


0/0 


+80%/-60% 




544 


1309 


5170 


12093 


w+w- 


38 


4 


16 


2 


zz 


8 


7 


36 


32 


Wev 


208 


60 


7416 


2198 


eeZ 


2 


2 


< 7 


< 6 


qq,q^t 


42 


45 


15 


17 


ti 








7 


7 


2-photon 


53 


53 


12 


12 


total background 


351 


171 


7509 


2274 


S/B 


1.5 


7.6 


0.7 


5.3 



Table 4: Expected numbers of events remaining at y / Sth = 260 GeV and ^/SpQ = 500 GeV, 
with unpolarized and with polarized beams, after sequential selection cuts have been applied. 
The entries in the form < N show the number of events corresponding to a single selected 
simulated event. 

evident from the table, the background can be greatly reduced for ^/s^, resulting in a 
very good signal-to-background ratio. At v /s^k, on the other hand, a large background 
from single- W boson production is left. For unpolarized beams, the resulting signal- 
to-background ratio is 0.7 While this would allow an unambiguous discovery of stop 
quarks (see Section |6.4|) , it is not a very good basis for precision measurements of 
the stop mass. Fortunately, the signal-to-background ratio can be greatly improved 
by using polarized beams. With an essentially right-handed electron beam and left- 
handed positron beam, the signal is enhanced, while most backgrounds are substantially 
suppressed. As a result, the signal-to-background ratio at — 500 GeV is improved 
from 0.7 to 5.3. 

We checked for possible supersymmetric backgrounds. The main concern is chargino 
pair production with the decay channel xt ~ * tib. We simulated a sample of these de- 
cays, consistent with our benchmark scenario, and found that the cuts listed in Table |3| 
completely eliminate this background source. 

With the results listed in Table [| for polarized beams, we can compute the observ- 
able Y and its statistical error, obtaining Y = 0.1082 ± 0.0034 with a relative error 
of 3.1%. The corresponding stop quark mass would be 

m h = (122.5 ± 0.19) GeV (3.5) 

where the uncertainty depends on the slope, dY/dM = —0.01755, at Y = 0.1082. 
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Without positron polarization, P(e + ) = 0, the precision of the measurement is reduced 
by roughly 20%, resulting in AY/Y = 3.7% and Am^ = 0.23 GeV. Even in this case 
the statistical error is rather small. 

It should be recalled that the production cross-section is a strong function of the 
mixing angle, so the statistical error Am^ will also depend on it. In our reference 
scenario, the light stop eigenstate is almost completely composed of the partner of the 
right-handed stop, t\ ~ tn, with the mixing angle cosOf = 0.01. While this scenario 
is preferred by electroweak precision data and the explanation of baryogenesis, an 
experimental analysis should consider all possible values for the stop mixing angle. For 
other values of cos^, the production cross-section can change drastically, depending 
on the beam polarization. As concrete example, we consider two larger values of cos 6f. 

cos %= 0.6: a L , 260 = 52 fb, a^ 500 = 194 fb, (3.6) 

ctr,260 = 39 fb, a Ri50 o = 148 fb, 

cos6> t ~ = 1.0 : er L ,260 = 169 fb, a L ,5oo = 577 fb, (3.7) 

0-R,26o = 6.9 fb, a Ri 5oo = 30 fb. 

Here cl/r^ stands for the stop production cross-section at center-of-mass energy E GeV, 
and with beam polarization combinations P(e~) = — 80%/P(e + ) = +60% and P{e~) = 
+80%/P(e + ) = —60%, respectively. If the stop is dominantly left-chiral, with | cos^| > 
0.5, the production cross-section is substantially larger for left-handed electron and 
right-handed positron polarization, opposite to the situation for a right-chiral stop. 
Therefore, for large values of | cos^l, it is better to use the beam polarizations P(e~) = 
— 80%/P(e + ) = +60%, even though one has to deal with much larger Standard Model 
background. The largest background, e + e~ — > Weu, amounts to about 12800 events 
at y/s = 500 GeV and £ = 200 fb - for this polarization. Nevertheless, due to large 
signal cross- sect ions, the resulting statistical error is still small, as summarized in Ta- 
ble [5], which demonstrates that, for all values of the stop mixing angle, one can measure 
the stop mass with a statistical error better than 0.3 GeV using our method and an 
appropriate choice of beam polarization. 

3.3 Iterative Discriminant Analysis 

A traditional, sequential-cut analysis was presented in the previous section. Often, 
more advanced multi-variate techniques can boost the sensitivity of a search. We 
investigated the efficacy of an Iterative Discriminant Analysis (IDA) for the purposes 
of measuring the stop quark mass based on the observable Y. 

The IDA method [|J is a modified Fisher Discriminant Analysis, the two main 
differences are the introduction of a non-linear discriminant function and iterations in 
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cos 9f 


P(e- 


) = -80%/P(e+ 


) = +60% 


P(e-) = +8 0%/P(e+) = -60% 


0.0 




0.69 




0.19 


0.6 




0.29 




0.28 


1.0 




0.14 




0.94 



Table 5: Statistical uncertainties Amj in GeV, for selected values of cos 9f and two opposite 
sets of beam polarization. The bold numbers indicate the best choice of beam polarization 
for the given value of the stop mixing angle. 

order to enhance the separation of signal and background. Two IDA steps have been 
performed. In order to have two independent samples for the derivation of the IDA 
function and for the expected performance, the signal and background samples were 
divided into two equally-sized samples. For this analysis, the same kinematic variables 
and simulated event samples as in the cut-based analysis are used, including the charm 
tagging flags F { . Before the mult i- variable analysis is performed, cuts on the input 
variables, so-call "tail cuts," are applied in order to improve the IDA performance. This 
is achieved by reducing the number of input events, and thus reducing the computa- 
tional time. From the distributions of the input variables for the signal and background 
events, the IDA method calculates a separating surface in the multi-dimensional pa- 
rameter space between signal and background events. The IDA output variable has a 
different shape for signal and background events, and therefore a cut on this variables 
is used to separate signal and background. In the first IDA step a cut is placed on this 
IDA output variable such that 99.5% of the signal efficiency are kept. The number of 
background events is largely reduced. From the smaller background sample and the 
99.5% remaining signal events again a new IDA output variable is constructed. The 
cut on the IDA output variable in this second IDA step defines the signal efficiency and 
the corresponding number of background events. Different working points are possible: 
they are defined by choosing a certain signal efficiency and obtaining the corresponding 
number of background events. The working point was determined by the expected error 
on mi x . The results of the IDA method with stop fragmentation are shown in Fig. |8| 
and Table |^ expressed as number of expected background events for each contributing 
background process. 

As before, in the channels where no event is left after the signal selection, an upper 
limit corresponding to one simulated event is given in the table. 

The IDA method achieves a significantly more powerful discrimination between 
signal and background than the analysis with conventional cuts. When allowing similar 
background levels as for the cut-based analysis in Table [|, signal efficiencies of e t h = 
0.387 for y/stt = 260 GeV and e p k = 0.416 for v /s^k = 500 GeV are obtained. 
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(a) yfi = 260 GeV (b) yfe = 500 GeV 




Efficiency Efficiency 



Figure 8: Performance of the Iterative Discriminant Analysis (IDA) for y/s = 260 GeV and 
y/s = 500 GeV. The plots show the remaining background event numbers for unpolarized 
beams as a function of the signal efficiency. 

With the resulting event numbers given in Table[| for P(e~) / P(e + ) = +80%/— 60%, 
the ratio quantity in eq. ( |2.3|) amounts to Y — 0.0648 ±0.0018 with a relative statistical 
error of 2.7%, translating into 

m h = (122.5 ± 0.17) GeV (3.8) 

where the uncertainty on the mass depends on the slope dY/dM = —0.01052. The 
higher signal efficiency and lower background achieved by the two-step IDA results in 
a slightly smaller statistical uncertainty (cf. Eq. fl3.5|) ). 



4. Experimental Systematics 

The high signal efficiency and low backgrounds achieved in both the cut-based analy- 
sis (Section [3~2"D and the IDA (Section deliver an excellent statistical precision - 
Amj < 0.2 GeV. It remains to investigate systematic uncertainties, which were the 
dominant contribution in the previous analysis of Ref . . We considered the following 
important sources of systematic errors: 

• detector calibration (energy scale) 
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V^th = 


= 260 GeV 


V^Pk 


- 500 GeV 




c = 


50 ftV 1 


C = 


= 200 ftv 1 


P(e-)/P(e+) 


0/0 + 


S0%/-60% 


0/0 


+80%/-60% 


Mi 


619 


1489 


9815 


22958 


W + W 


11 


1 


< 8 


< 1 


ZZ 


< 2 


< 2 


20 


18 


Wev 


68 


20 


1719 


510 


eeZ 


3 


2 


< 7 


< 6 


qq,q^t 


16 


17 


18 


21 


tt 








1 


1 


2-photon 


27 


27 


294 


294 


total background 


127 


69 


2067 


851 


S/B 


4.9 


22 


4.7 


27 



Table 6: As in Table ||, expected numbers of events remaining at ^/s t h = 260 GeV and 
y/Spk = 500 GeV, with unpolarized and with polarized beams, after the IDA has been applied. 



• charm tagging 

• hadronization / fragmentation 

• neutralino mass 

• luminosity measurement 

• beam energy spectrum 

• background estimate 

The first four sources pertain to the signal efficiency We discuss these sources in detail 



in the context of the sequential-cut analysis detailed in Sec. 3.2 first, and then briefly 



report the results obtained from the IDA method described in Sec. 
4.1 Systematics for the Sequential-Cut Analysis 

Many of the kinematic quantities used in these selections depend on a correct calibration 



of the calorimetry. Based on experience from LEP [ 27| , we assume an uncertainty of 1% 



on the overall energy scale, which is rather pessimistic for a future ILC detector. We 
scaled simultaneously all kinematic quantities through a range of ±6% and observed 
correlated shifts in the overall selection efficiency at the two center-of-mass energies. 
In particular, the px cut is sensitive to this kind of scale error, prompting us to tune 
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the cut at Ws p k to achieve a minimal residual uncertainty for the ratio of efficiencies, 



as discussed in Section 3.2 



Figure [5] shows how the selection efficiencies change as a function of the scale factor. 
Using our optimized px cut shown in Fig. |5|, one sees a parallel behavior at threshold 
(upper solid line) and at peak (lower solid line). This leads to a very good cancellation 
for the ratio of efficiencies, as shown by the solid line in Fig. |1(]. If we had optimized 
for efficiency only, then we would have used nearly the same px cuts at the peak as 
we use at threshold. However, this would have given a rather different dependence 
on the scale, as indicated by the dashed line in Fig. |9|, and therefore a much stronger 
dependence of the ratio of efficiencies on the scale, as shown by the dashed line in 
Fig. [l^. With our best cuts, an uncertainty of ±1% on the calorimeter energy scale 
translates into an uncertainty of less than 0.6% on the ratio of efficiencies. 
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Figure 9: Variation of the selection efficiencies (eth and e p k) with an overall energy scale 
factor. The two solid lines show the variation obtained with our nominal cuts, at y^Sth 
and ^/Spk. The dashed line shows what we would obtain if we applied a looser px cut at WSp^. 



The efficiency for track reconstruction should be very high at an ILC detector. 
However, there is always an uncertainty in the value for that efficiency, which we took 
to be about 0.5%. We propagated this uncertainty to the cut on iV" tracks , since a random 
loss of tracks changes the shape of the distribution of Arracks- Since our cut is quite 
loose, however, the resulting uncertainty on the ratio of efficiencies is negligible. 

(c) 

Knowledge of the efficiency for charm jets for a given cut on F£ is not easy to 
obtain. Based on the work described in Ref. p6 |, we assumed an uncertainty of 0.5% 
on the charm efficiency. Although one might expect this uncertainty to be correlated 
between the two center-of-mass energies, we assumed no strong correlation and assign 
an uncertainty of 0.5% for the ratio of efficiencies. 
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Figure 10: Relative variation of the ratio of efficiencies (e p k/eth) with an overall energy 
scale factor. The solid line shows a very small variation, given our nominal cuts on pr, to be 
compared to a much larger variation if we had used looser pt cuts meant to maximize the 
efficiency 



The observable Y depends on the integrated luminosity at both center-of-mass en- 
ergies. Traditionally, the luminosity is measured using Bhabha scattering, for which 
highly accurate theoretical cross-sections are available. The limiting systematic uncer- 
tainty for the LEP detectors comes from the acceptance of the luminosity calorimeters. 
Such an uncertainty would essentially cancel in the ratio of luminosities. Alternatively, 
one could define an effective luminosity through another clean QED process, such as 
e + e~ — > /U + /i~, for which there is essentially no theoretical or experimental systematic 
effect. The precision of the ratio of luminosities would come from the statistical uncer- 
tainty on the number of events recorded, which we estimate to be about 0.4%; 
this is the figure we use in this study. 

Apart from instrumental issues such as the energy scale, track reconstruction ef- 
ficiency, charm tagging efficiency and the measurement of the integrated luminosity, 
the estimate of the signal efficiency will depend on the modeling of the signal itself. 
While the production of a pair of stop quarks is well understood and can be modeled 
accurately, the non-perturbative aspects of the formation of stop hadrons which then 
decay into two or more jets are more problematic 1 . We have attempted to account 
for this fundamental difficulty by varying the parameter which controls the fragmen- 
tation function in our simulations. We used PYTHIA and the Peterson fragmentation 
function, with values of the fragmentation parameter reported by the OPAL Collabo- 
To be specific, we took e c = -0.031 ± 0.011 and e b = -0.0050 ± 0.0015, 

2 



ration p8[ . To be specific, we took e c = —0.031 ± 0.011 and e b = 
and propagated e b according to the assumption that ej = e b (m b /mfj 



We varied e c and ej independently, and measured the efficiencies at the two center- 



1 Earlier analyses such as Ref. |jj neglected this important problem. 
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of-mass energies. The impact of varying e c is small. For variations of the stop quark 
fragmentation, however, we find that the variation of the efficiencies is rather different, 
so the desired cancellation of systematic uncertainties is not achieved. In fact, most of 
the systematic uncertainty comes from the cut on px, and has an opposite sense at v /s^h 
and y/~Spk. 

The range in e b used in our simulations is quite broad. The more advanced measure- 
ments of 6-quark fragmentation from ALEPH [30] and OPAL [31] give more constrained 
values: e b = -0.0031 ± 0.0006 (ALEPH) and e b = -0.0041 ± 0.0004 (OPAL), using 
rather different methodologies. On the basis of these measurements, one could argue 
that our assumed variation in e b is too large by a factor of two. 

Rather than relying on LEP measurements to predict stop quark fragmentation, 
we investigated the potential of ILC data to constrain the fragmentation. We already 
noted that most of the sensitivity to stop quark fragmentation comes from the cut 
on pt'-, however, the change in the shape of the px distribution is small. (The fact that 
the quantities chosen for cuts are insensitive to e t ~ is a strong point of the analysis, of 
course.) We examined other kinematic quantities and found a few which exhibit clear 
changes in shape when we vary ej. Four examples are shown in Fig. |Tl|. The M v \ s / \fs 
distribution shows pronounced shifts as a function of e t ~. Given an accumulation of a 
few xlO 4 events at ^/s^k, one can show that the mean of this distribution alone would 
allow a differentiation of our three values e b = —0.0050 ±0.0015 at more than ten sigma 
(statistical uncertainty only). If the energy scale uncertainty were a problem, then one 
could normalize M v i S to E V - 1S — a clear distinction between the three distributions is 
visible near the peak of M v - ls / E w - ls . The energy of the third jet, when it exists, shows a 
good sensitivity to e 4 -. (Recall that the jets are ordered in decreasing energy.) Better, 
perhaps, is the smaller of the two di-jet invariant masses formed by combining this third 
jet with the first and second jets. Although these considerations are not equivalent to 
a full study of a possible measurement of the stop fragmentation, they do indicate 
that a good measurement should be possible, well beyond the extrapolation of LEP 
results on e b to e t ~ and all the attendant assumptions behind such an extrapolation. 
On this basis, we judge that the uncertainty on the stop fragmentation would be no 
larger than one-fourth of the uncertainty obtained by comparing simulations with e b = 
-0.0035, -0.0050 and -0.0065, which corresponds to Ae f = (2.5 x 10" 6 ) /4 = 0.6 x 10~ 6 . 

Another empirical quantity which induces an uncertainty on the selection efficiency 
is the mass of the neutralino, rn^o. The mass difference — m^o clearly impacts the 
kinematic distributions, so the efficiency estimated from the simulation depends directly 
- and dramatically — on m^o. We simulated a sample with m^o = 108.2 GeV, which is 
one GeV higher than our default value. The relative change in the selection efficiencies 
is roughly 10%. Since the changes are parallel, the ratio of efficiencies change by 
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Figure 11: Changes in kinematic distributions at ^/s p k for different assumed values of ej, 
which are related directly to the listed values of e& through = e& (mb/m^) 2 . The solid points 
with error bars show the distribution with e& = —0.0050, our default choice. The thick, 
unshaded red histogram shows e& = —0.0035, and the thin, shaded blue histogram shows 
e b = -0.0065. 



only 2.8%, once again illustrating the robustness of this method. Other studies have 
shown [|J that m^o can be measured with an accuracy of 0.3 GeV or better, so we 
assign an uncertainty of 0.8% due to the unknown neutralino mass. 
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The predicted cross-sections depend on the beam energy and the beam energy 
spectrum. Due to beamstrahlung and other effects, the mean energy can be significantly 
lower than the peak value. While we used CIRCE for taking this fact into account, 
the question remains how well a program such as CIRCE can be validated using real 
data. This question has been addressed by several authors, using, for example, Bhabha 



scattering and radiative returns to the Z pole p2" |. The studies indicate that models 
for the spectrum and the beam energy can be constrained directly from the data to an 
accuracy on the order of 0.1 GeV. We include this uncertainty as a direct uncertainty 
on m^, but not on the observable Y . 

Finally, we must consider uncertainties on the estimated contributions from back- 
ground processes. The SM backgrounds fall naturally in two categories: two-photon 
interactions, which are difficult to predict, and the others, which involve high-p^ electro- 
weak processes, for which direct theoretical calculations are reliable. We also consider 
background contributions from the production of other supersymmetric particles. 

Two-photon interactions cannot be fully described by perturbative QCD, and so 



phenomenological models are required [33]. These must be tuned to match real data, 



which is difficult due to the fact that most two-photon scattering events emit particles 
that are lost down the beam pipe. Parameters pertaining to the softest interactions are 
difficult to constrain; fortunately, such interactions are easily eliminated by our cuts 
on p T , iV tracks and E vis . Many of the events coming in at higher p T can be described using 
models with a basis in perturbative QCD. The investigations of the photon structure 
functions by the LEP Collaborations illustrate the procedure of tuning parameters 
and confronting the models with real data, leading to interesting conclusions about 
the success of the various models |3j§, 1351 BBf. It is not straightforward to translate 



those conclusions into constraints on our two-photon background, although Figs. 19, 21 
and 23 in the report from the LEP Working Group |$6| and Figs. 5-7 in the OPAL 
paper |55 are quite relevant for our study, and indicate that modeling the tails of the pt 
distribution at the 20% level should be possible. Assuming that the study of two-photon 
interactions would be greatly extended at the ILC, we assign a 20% uncertainty to the 
background estimate for two-photon interactions. The resulting relative uncertainty on 
the Y observable is 0.8%. 

The dominant background is e + e~ — > Weu, according to Table |] (and Table ^[). 
A precise prediction of this background requires accurate measurements of this process 
combined with the calculation of higher-order radiative corrections. While a complete 
NLO calculation of that process is missing, a recent result for the related process of 
W pair production |37|] suggests that a NLO calculation of Wev is feasible within the 
next years with an error remaining well below 0.5%. The impact on Y is negligible, on 
the order of 0.1%, relative. 
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error on 


relative shift in expected signal yield (%) 


error on 


variable 


variable 


■s/s^ = 260 GeV 


^s~^ = 500 GeV 


Y {%) 


energy scale 


1% 


3.7 


3.1 


0.6 


tracking efficiency 


0.5% 




negligible 




charm tagging efficiency 


0.5% 




taken to be 0.5 




luminosity 




0.4 


0.2 


0.4 


charm fragmentation 


0.011 


0.3 


0.8 


0.6 


stop fragmentation 


0.6 x lO" 6 


0.6 


0.2 


0.7 


neutralino mass 


0.3 GeV 


3.8 


3.0 


0.8 


background estimate 




0.8 


0.1 


0.8 



Table 7: Evaluation of experimental uncertainties on Y, for the sequential-cut analysis. 
The last column gives the relative uncertainty on Y. 



A summary of the experimental systematic uncertainties for the sequential-cut anal- 
ysis is shown in Table 0. A good cancellation of experimental systematics is obtained, 
except for the stop quark fragmentation uncertainty and the background estimation. 
The goal of the new method is therefore fairly well achieved with this set of sequen- 
tial cuts. The implications for the measurement of the observable Y and the inferred 
mass mfj are discussed in Sec. ||. 

4.2 Systematics for the Iterative Discriminant Analysis 

We evaluated the impact of the sources of systematics listed on page [H] in a manner 
similar to the methods of Sec. [Q[ We scaled all kinematic inputs to the IDA ac- 
cording to an overall energy scale uncertainty. The systematic uncertainty from the 
number of tracks is assumed to be negligible. The variations in the charm and stop 
quark fragmentation functions were evaluated as before. The sensitivity to m^o and 
the uncertainty on the background estimate were evaluated precisely as above. The 
luminosity uncertainty is, of course, the same as in the sequential-cut analysis. 

The resulting systematic uncertainties are listed in Table We observe a much 
larger uncertainty coming from the scale uncertainty as compared to the sequential-cut 
analysis (cf. Table [7|). With multi-variate methods such as the IDA, it is difficult to as- 
certain what role any given quantity plays in the final output variable, so no dissection 
of the IDA to reveal the sensitivities to the energy scale is possible. Furthermore, one 
cannot tune the operation of the IDA in order to balance efficiencies for each quantity, 
as we did for thrust T and event-p^ in the sequential-cut analysis. For this kind of 
precision measurement, it would appear that the better discrimination of signal and 
background provided by the IDA is of limited value in light of the larger and uncontrol- 
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error on 


relative 


shift in expected signal yield (%) 


error on 


variable 


variable 


y/sth = 


260 GeV 


= 500 GeV 


Y (%) 


energy scale 


1% 


3.4 




1.3 


2.3 


tracking efficiency 


0.5% 






negligible 




charm tagging efficiency 


0.5% 






taken to be 0.5 




luminosity 




0.4 




0.2 


0.4 


charm fragmentation 


0.011 


0.1 




0.6 


0.5 


stop fragmentation 


0.6 x lO" 6 


0.1 




0.8 


0.7 


neutralino mass 


0.3 GeV 


3.7 




1.6 


2.2 


background estimate 




0.3 




0.2 


0.1 



Table 8: Evaluation of experimental uncertainties on Y, for the IDA. The last column gives 
the relative uncertainty on Y. 



lable sensitivity to experimental sources of systematic uncertainty. Nonetheless, when 
performing a measurement with real data, one would welcome an alternative analysis 
in order to check the robustness and stability of the measurement. 

In Section B71, we show the power of the IDA in the discovery of a light stop quark. 



5. Theoretical Uncertainties 

The inference of the stop mass from stop cross-section measurements requires precise 
theoretical calculations for the cross-sections. The stop production cross-section re- 
ceives large corrections in particular from QCD gluon exchange between the final state 
stops. Near threshold, when the stop quarks are slowly moving, these effects become 
very large, which is the well-known Coulomb correction The NLO QCD correc- 
tions to stop production have been computed several years ago |15| and it was found 
that the corrections range between about 10% at high energies and up to 100% near 
threshold. This shows that higher-order corrections are crucial. 

Over the last few years, sophisticated techniques have been developed for calculat- 
ing top-pair production at NNLO |39|]. Near threshold, they include resummation of 



terms of order 0(a s /v) for the low velocity v of the top-quarks. For the production 
of scalar quarks, similar calculations are not yet done. However, one can use partial 
results to obtain a reliable estimate of the uncertainty of the NNLO corrections. At 
NNLO order, several contributions enter in the computation. The largest effect near 
threshold arises from the Coulomb correction. The impact of the Coulomb corrections is 
calculated through NNLO order [40|] , including resummation via non-relativistic QCD. 



-27- 



Technically, here the non-relativistic Schrodinger equation is used for computing the 



Coulomb effects 41 



Similar to the case of top pair production, it is found that the NNLO term to stop 
pair production is of similar order of magnitude as the NLO term, i.e., the perturbation 
series is converging rather slowly. From the behavior of the perturbation series and 
the remaining scale dependence, the size of the missing higher-order contributions is 
estimated to be around 7% at 260 GeV and 2.5% at 500 GeV. 

However, we want to point out that several improvements to this straightforward 
approach could be made. Besides the large Coulomb- type corrections of order 0(a s /v), 
there are also potentially large logarithmic contributions 0(log(a s /t> )). They can be 
resummed with more sophisticated techniques, for instance velocity non-relativistic 
QCD [Q. Using the results of Ref. ||3| for the NLO corrections to squark pair produc- 
tion, it is found that the uncertainty with respect to the NLO computation is reduced 
significantly. A similar improvement can be expected at the NNLO level. In addition, 
instead of directly computing the total cross-section near threshold, one can describe 
it through moments \M\ that avoid the non-perturbative contribution of stopponium 
bound states that can form just below the nominal stop-pair threshold. With these 
refinements it is expected that the theoretical uncertainty can be brought down by 
a factor of two (however the actual calculation remains for the future). So here an 
uncertainty of 3.5% at 260 GeV and 1% at 500 GeV are assumed. 

Besides the QCD corrections, the electroweak corrections need to be considered. 
The NLO electroweak corrections have been computed and found to amount to 
several per-cent. While they need to be taken into account, the NNLO corrections are 
expected to be much less than 1%, with the exception of leading initial- and final-state 
QED corrections that can easily be resummed to higher orders. 

Combining the two errors, a total theoretical error of 4% at 260 GeV and 1.5% at 
500 GeV can be assigned. Pessimistically, we add these two uncertainties linearly, and 
assign a theoretical uncertainty of 5.5% for the quantity Y. 



6. Results and Implications 

We derive the expected measurement error on the stop quark mass and discuss the 
implications for particle physics predictions of the relic density of dark matter. We also 
discuss the luminosity needed to discover a light stop quark in this scenario, using the 
IDA method. 

6.1 Precision on the Stop Quark Mass 

A final assessment of the achievable precision on the stop mass will be based on the 
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statistical and all systematic uncertainties. Table ^| summarizes these uncertainties for 
the observable Y defined by Eq. ( |2.3|) . One sees that the IDA method achieves a smaller 
statistical uncertainty on Y at the cost of a larger experimental systematic uncertainty. 
It would be important, in a measurement with real data, to implement two methods 
as we have done here, and check the consistency of the results. 

The stop quark mass is inferred from the measured values of the observable Y 
following the example described in Sec. |||. The differing efficiencies for the sequential- 
cut and IDA methods lead to different central values for Y and for the slope dY/dM 
at that point. The inferred uncertainties on the stop quark mass are summarized in 



Table |10| and are similar for the two analyses. We conclude that the stop quark mass 
could be measured with an uncertainty of Am^ = 0.42 GeV. 

We investigated the dependence of the measurement error on the integrated lu- 
minosity. There is very little change in the statistical uncertainty if we increase the 
luminosity on peak, but the variation with the luminosity at y/s^ = 260 GeV is in- 
teresting - see Fig. |12| The experimental uncertainty is dominated by the statistical 
contribution, so a decrease in the luminosity from our assumed value of £ t h = 50 fb -1 
has a significant impact. On the other hand, the theoretical uncertainty is very large 
by comparison, so increasing £ t h hardly improves the total error on . A luminosity 
in the range 30 fb" 1 < C th < 80 ftT 1 would appear to be optimal, for this analysis. 

The dominant uncertainty comes from the theoretical calculation of the signal 
cross-section. As discussed in Sec. this uncertainty comes mainly from higher-order 
corrections which are not easily summed at threshold. The estimate of this theory error 
relies on present computational techniques and some expectations on how they might 
improve in the future. However, the progress in calculations of radiative corrections 
can not really be predicted, so the assumed value for the theoretical uncertainty at the 
time when ILC is running might well be somewhat different than the value reported in 
Table |. In particular, history has shown that people working on loop computations 
often overcame big problems with unexpected ingenuity, in order to be able to make 
most of precise measurements. Therefore, in the following, the combined error in 
Table ^ will be taken as a conservative estimate. If one were to set aside the theoretical 
error on the cross section, then the total experimental error is quite small, amounting 
to 3.5%-4.3% on Y. In this case, the error on the stop quark mass would be a little 
larger than Am^ = 0.2 GeV. 

6.2 Comparison with Previous Results 

A previous study investigated the potential of the ILC running at y/s = 500 GeV to dis- 
cover a light stop quark and measure its parameters JTJ. It was assumed that 250 fb -1 
would be taken at two beam polarization combinations: P(e~)/P(e + ) = +80%/ — 60% 
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error source for Y 


sequential cuts 


IDA method 


detector effects 


0.9 


2.4 


charm fragmentation 


0.6 


0.5 


stop fragmentation 


0.7 


0.7 


neutralino mass 


0.8 


2.2 


background contribution 


0.8 


0.1 


sum of experimental systematics 


1.7 


3.4 


statistical 


3.1 


2.7 


sum of experimental errors 


3.5 


4.3 


theory for signal cross-section 


5.5 


5.5 


total error AY 


6.5 


7.0 



Table 9: Summary of relative statistical and systematic uncertainties (in percent) on the 
observable Y. 





measurement error Am^ (GeV) 


error category 


sequential cuts 


IDA method 


statistical 


0.19 


0.17 


sum of experimental systematics on Y 


0.10 


0.21 


beam spectrum and calibration 


0.1 


0.1 


sum of experimental errors 


0.24 


0.28 


sum of all experimental and theoretical errors 


0.42 


0.44 



Table 10: Estimated measurement errors (in GeV) on the stop quark mass 



and —80%/ + 60%. Measurements of the stop squark production cross sections at these 
two polarizations are sufficient to deduce the mixing angle and mass of the stop squark. 
A host of systematic uncertainties was considered, with the conclusions that the absolute 
cross-sections could be measured to 1.3%-2.4%, dominated by experimental systematic 
uncertainties (the statistical uncertainty was 0.8%). Under the given theoretical sce- 
nario, the total error on the stop quark mass was estimated 2 to be Am,^ =1.2 GeV. 

Our theoretical scenario coincides with the one studied in Ref . [[[[ , and the method 
proposed here leads to a total error on the stop quark mass that is more than two 
times smaller: Am^ = 0.42 GeV, even though a much smaller integrated luminosity is 
assumed. This improvement is certainly valuable, and is quite helpful for the calculation 
of the relic density, as we discuss below. We would like to point out, however, that the 

2 Note that the error of 1.2 GeV is slightly larger than reported in Eq. (17) in Ref. [Q, since we are 
using the scenario from Ref. || with large slepton masses, leading to a larger neutralino mass error, 
which in turn increases the stop mass uncertainty. 
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luminosity at threshold (fb~ ) 



Figure 12: Decrease of the statistical uncertainty (blue line), total experimental uncer- 
tainty (red line) and total uncertainty on mr (black line), as a function of the integrated 
luminosity £ t h at y / Sth. 

basis for the experimental analysis in Ref. [|lj] differs significantly from what was used 
for the present analysis. In particular, the fragmentation of the stop squark and of the 
charm quark produced in its decay was not simulated in Ref. |l| , leading to very different 
signal characteristics which are not realistic. For example, the number of jets was almost 
always two, which contrasts starkly with the present study in which typically one or two 
additional jets are found due to the process of fragmentation. Thus the requirement 
in Ref. |TJ of exactly two jets, which is very effective at suppressing the simulated 
Wev background, leads in practice to a very low efficiency and large sensitivity to the 
modeling of the fragmentation practice. Ref. |lj assumed a 1% systematic uncertainty 
on the absolute efficiency due to fragmentation, which dominates the total uncertainty 
on the cross-section measurement. As we have seen in the present study, this is likely to 
be significantly underestimated. Furthermore, no serious assessment of the theoretical 
uncertainty on the cross-section was given in Ref. [[!]; this uncertainty was assumed to 
be negligible. As we discuss in Section [| above, this assumption is not justified, and 
with existing techniques we estimate the theoretical uncertainty to be as large as 2.5% 
at a/s = 500 GeV, which we assume will be reduced to ~ 1% by the time the ILC is 
running. 

We evaluated the efficiency of the selection in Ref. |I[ using the signal samples 
generated for the present study. We find that the requirement A7j et s = 2 delivers an 
efficiency of only 7.3%, to be compared to ~ 18% reported in Q. This would increase 
the statistical uncertainty on the cross-section measurement from 0.8% to 1.2%. If we 
relax the cut to iVj Cts > 2, the efficiency climbs to 21%, though clearly the background 
would become too large with this cut. We evaluated the scale dependence and found it 
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to be negligible, since there are relatively few energy-based cuts in IJ. We evaluated the 
systematics for the stop and charm fragmentation, and found a very large sensitivity 
to the fragmentation, on the order of 5%, due to the requirement that iVj Cts = 2 only. 

We conclude that the uncertainty on the stop quark mass, Am^ = 1.2 GeV, re- 
ported in Ref. was underestimated, so that our present result Aro^ = 0.42 GeV 
represents a major step forward. 

6.3 Implications for Relic Density Calculation 

Precise measurements of supersymmetric particle properties at the LHC and ILC can 
be used to compute the dark matter relic abundance so as to compare with cosmological 
observations. If stop-neutralino co- annihilation is relevant, as in the scenario studied 
here ||, it is important to measure the stop-neutralino mass difference very precisely. 
The extraction of the neutralino properties, in particular the lightest neutralino mass, 
is studied in detail in Ref. ||. It is found that a high precision of Am^o w 0.3 GeV for 
the lightest neutralino mass can be achieved at the ILC, and also the other neutralino 
parameter can be inferred rather well. 

The limiting factor in the accuracy of the dark matter estimation is therefore the 
precision of the measurement of the scalar top quark mass. As discussed in the previous 
section, an older study using cross-section measurements at y/s = 500 GeV found 
Anifj = 1.2 GeV and | cos 6^1 < 0.077 for the stop mass and mixing angle [I], ||. Based 
on these expected experimental results, the relic dark matter density is computed with 



the codes described in Ref. |TT|, B6f. Fig. 13 shows the result of a scan over the MSSM 



parameter space. The scattered gray dots indicate the region allowed by the collider 
experimental uncertainty, as a function of the measured stop mass. The horizontal 
bands depict the relic density as measured by WMAP || with one and two standard 
deviation errors. Here, ^cdm is the ratio of the dark matter energy density to the critical 
density p c = 2Hq/(8itG^), with the Hubble constant H = h x 100 km/s/Mpc and 
Newton's constant G^. At the la level, the astrophysical observations lead to 0.103 < 
^cdm^ 2 < 0.116. With a stop mass measurement error of = (122.5 ± 1.2) GeV, 
the relic density can be predicted to 0.082 < Ocdm^ 2 < 0.139 at the la level. With 
the new result of this work, Am^ = 0.42 GeV, the relic density can be computed 
much more precisely, yielding the result 0.096 < ^cdm^ 2 < 0.124. This precision is 
very comparable to the direct WMAP measurement, as indicated by the black dots in 
Fig. 0. 

As pointed out above, the estimate of a stop mass error of Am^ = 0.42 GeV is 
based on a rather conservative evaluation of systematic errors. In particular, this value 
is dominated by the conjectured theory error on the prediction of signal and background 
cross-sections. If on the other hand, with progress in calculation methods, the theory 
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Figure 13: Computation of dark matter relic abundance ^cdm^ 2 taking into account esti- 
mated experimental errors for stop, chargino, neutralino and Higgs sector measurements at 
future colliders. The dots correspond to a scan over the la (Ax 2 < 1) region allowed by the 
experimental errors, as a function of the measured stop mass, for a mass measurement error 
of 1.2 GeV (light gray dots), 0.42 GeV (dark gray dots) and 0.24 GeV (black dots). The 
underlying scenario used as input is indicated by the red (light) star. The horizontal shaded 
bands show the la and 2a constraints on the relic density measured by WMAP. 

error could be reduced to a sub-dominant level, the remaining statistical and systematic 
experimental errors would give a stop mass error of Amf = 0.24 GeV for the cut-based 
analysis and Am^ = 0.27 GeV for the IDA. The amelioration of the prediction for the 
dark matter relic density due to this improvement in stop mass precision is illustrated 
in Fig. [L3|. 

For this accuracy of the stop mass measurement, the uncertainty of the dark matter 
prediction becomes limited due to the expected experimental errors in the lightest 
neutralino mass and mixing angles, which we have taken from Ref. 0. As a result, 
taking an error of Am^ = 0.24 GeV for the stop mass, we find 0.099 < f2cDM^ 2 < 
0.121, which is only a small improvement in the precision of the dark matter density 
prediction with respect to Am^ = 0.42 GeV. 

6.4 Discovery of the Light Stop Quark 

The main focus of this paper is the measurement of the stop quark mass. It is inter- 
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esting, nonetheless, to consider the utility of these selections for discovering the light 
stop quark at the ILC. The IDA-based selection, in particular, achieves a very low 
background and a high efficiency - see Table | in Subsection |3T3"| . 

We examined this issue assuming that the ILC collects data at ^Spk = 500 GeV, 
with unpolarized beams, as one might expect at start-up. The signal cross-section 
for this scenario is cr f - = 118 fb. The nominal IDA selection efficiency is e = 0.416 
and the background cross-section for unpolarized beams is o\, = 10.3 fb. Tightening 
the selection to reduce the background improves the sensitivity of the analysis only 
very slightly. This information allows a calculation of the expected tail probability 
or p-value as a function of integrated luminosity, C Specifically, we computed the 
p-value setting the hypothetical number of observed events equal to the mean of the 
corresponding Poisson distributions (signal and background), as a function of C. The 
result is shown in Fig. [14] by the thick red line. The black dots on the line show 
hypothetical integral numbers of observed events, starting at N — 1 for C = 0.02 fb -1 . 
The plot clearly indicates that a luminosity of only C ~ 240 pb _1 would produce eleven 
observed events, on average, and the significance of ten signal events over the expected 
background would be more than 5cr. The uncertainty on the background estimate and 
the signal efficiency have a negligible impact on this result. 

This example applies only to our given scenario, with = 122.5 GeV, m^o = 
107.2 GeV and cos#£ = 0.01. Further investigations would be needed in order to 
understand how well this IDA selection would perform for other mass and mixing 
combinations. 




Figure 14: p- values as a function of integrated luminosity C. The black dots on the line 
show hypothetical integral numbers of observed events, starting at N = 1 for C = 0.02 fb -1 . 



7. Summary 



A new method for a precise measurement of the stop quark mass has been described, 
based on the ratio of yields at the peak stop quark pair production cross section, and 
near threshold. This ratio is far less sensitive to experimental uncertainties than other 
methods, leading to a very low estimated uncertainty, still dominated by the statistical 
uncertainty and the theoretical uncertainty (which is also present for other methods 
based on a cross-section measurement). We studied a specific scenario in detail, with 
an emphasis on analysis techniques and systematic uncertainties. We placed special 
emphasis on the modeling of the stop quark and charm fragmentation uncertainties, 
and suggest how fragmentation models could be constrained with data taken at the 
ILC Previous studies had not considered this source of uncertainty. This method is 
general, and could be applied to other species, provided an accurate prediction for the 
excitation curve is or can be available. For weakly interacting particles, such as staus, 
the theoretical uncertainty is much smaller and the advantage of the new method would 
be even more impressive. 

An important part of our studies is the use of multi-variate methods to isolate a very 
clean stop quark signal. For this we utilize the Iterative Discriminant Analysis (IDA) 
used previously at LEP. It is interesting that a carefully-tuned set of sequential cuts 
achieves a much smaller systematic uncertainty, and hence a better overall result for 
the stop quark mass measurement in this method. The superior background rejection 
of the IDA, however, is extremely useful when searching for a stop signal, and we give 
an illustration for y/s = 500 GeV, which shows that a five-sigma significance could be 
obtained by the IDA selection with only 240 pb _1 . 

The reduction of the uncertainty on the stop mass from about Am^ = 1.2 GeV 
in Ref. |3|] to Am^ = 0.42 GeV in this analysis is crucial for testing theoretical 
explanations of the dark matter relic density in the light-stop co-annihilation scenario. 
With these new results, the theoretical calculation has an accuracy equal to the two- 
sigma uncertainty of the WMAP measurements. The remaining uncertainty is no longer 
dominated by Am^. 
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